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Abstract 

B K*TT and B — > pK decays involve the same quark- level processes as 
B ttK. Analyzing the measurements of the former decays might be able to 
shed additional light on the new-physics hints in the current B — > ttK data. 
We perform fits to — > K*it and B — > pK decays, and find that the data 
can be accommodated within the standard model. However, this agreement 
is due principally to the large errors in the data, particularly the CP-violating 
asymmetries. If the errors on the B K*tt and B — > pK observables can 
be reduced, one will have a clearer sense of whether new physics is present in 
these decays. 
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In recent years, experimental data m B ^ nK decays have shown some dis- 
crepancies with the predictions of the standard model (SM), hinting at new-physics 
(NP) contributions in the decay processes If]. This so-called tcK puzzle has received 
a great deal of attention. The NP presumably enters the short-distance pieces of the 
decay amplitudes. It should therefore affect other closely-related processes as well. 
In particular, the final-state VP counterparts of B irK decays {V is a vector 
meson, P is a pseudoscalar @) - B ^ K*7i and B — > pK decays - involve essentially 
the same quark-level processes. It is thus of interest to see whether the current 
experimental data of these decays also present similar hints of NP. (Even within the 
SM, these decay modes are useful in extracting the weak phase 7 [SI El E]-) 

The eight VP decay modes can be readily divided into two sets: (I) B — » fC*7r 
decays and (II) B — > pK decays, each of which involves a different set of flavor 
amplitudes from the other. As is usually done, in our analysis we neglect the 
dynamically-suppressed exchange and annihilation amplitudes. In addition, the 
quantum chromodynamics (QCD) penguin amplitude has three pieces, and can be 
written as 

p' = v:,v^,p: + v:,v,p^ + v*v,p; 

= v:,vup'u-p'c) + v:mH-p':) . (1) 

(The primes on the amplitudes indicate h ^ s transitions.) In writing the sec- 
ond line, we have used the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix. In the following, all diagrams (except P^') are redefined to absorb the mag- 
nitudes of the CKM matrix elements which multiply them, i.e., P^ = |^*^Vu<(|P^, 
PL = \V:,V^s\{PL - P'c). PL = \V:M{PI - P;), etc. Because « \V:,Vl 

it is expected that \PL\ ^ \Ptc\- To begin with, we keep six types of diagrams: 
the color-allowed tree, denoted by T'; the color-suppressed tree, C"; two QCD pen- 
guin amplitudes, P^ and P^^,; the color-allowed EW penguin, P'^w'^ ^md the color- 
suppressed EW penguin, P'^w W^e further use a subscript P or for each diagram 
to indicate which final-state meson contains the spectator quark of the B meson. 
Explicitly, we have for the decay modes in Set (I): 
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^In the latest update of the ttK puzzle, it was seen that, although NP was hinted at i? ttK 
decays, it could be argued that the SM can explain the data, see Ref. ,1^ 

■*There are a great many articles on i? ^ VP decays, though few specifically on B ^ K*tt and 
B pK . For general papers on B ^ VP decays, see Ref. [2]. For analyses more specific to the 
B K*iT and B pK decays, see Refs. [3] and [1]. 
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and for the decay modes in Set (II): 
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A~+ — p' _ p' pH _ T' p'^f -I P"^ 

^11 — ^tcy ^ucy^ ^v^ ^ ^-^Ewy ' 
\/2A55 = ~Ptcy + P'ucy^^"' + -^w.p ~ C'pC*''' + -PEwy ! (3) 

In the amplitudes A^i (a G {/, //} and j G {0, +}) appearing in Eqs. ([2]) and ([3]), 
the subscript a denotes the decay set and the superscripts ij denote the electric 
charges of the vector and pseudoscalar mesons, respectively. Here we have explicitly 
written the weak phase associated with each amplitude (including the minus sign of 
V^^yts in the QCD and EW penguins), while the diagrams contain the strong phases. 

In fact, there are relations among the various diagrams. The ratios of Wilson 
coefficients Cio/cg and C2/C1 are equal, to a good approximation (the difference is 
only about 3%). In the limit in which this equality is exact, the relations read [S] 

P'ewm = -\- [T',i + Pi^M') M'^M, 

Z Ci 
Z Ci 

Now, the diagram does not appear in the -B-decay amplitudes. Thus, in order to 
use these relations, an approximation must be made. We consider two possibilities. 
In the first, we neglect all factors of P^ and P^^. In the second, we replace the P^ 
in the above relations by P^^. Both approximations are justified because |P^| and 
\Puc\ are expected to be small. In either case, the relations relate some diagrams 
appearing in P — *■ K*^ decays to some in P — pK decays. That is, the fits 
must be performed to measurements of both decays simultaneously. We call the 
fits corresponding to the two approximations Fitl and Fit2. Also, throughout this 
paper, the strong phases of the amplitudes are defined with respect to Tp. 

Current experimental measurements of the observables in P — * K*7c and P — ^ 
pK decays are given in Table [H In the first approximation, where P^^ is neglected, 
the direct CP asymmetries (Acp^s) of K*^n^ and p^K^ are identically zero within 
this framework. Therefore, it is of no use including them in the functions as they 
impose no constraint on the theory parameters. We are then left with 15 and 17 
observables in Fitl and Fit2, respectively. 
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Set 




Mode 


BR (xlO-6) 


Acp 


Sep 


(I) 


B+ 




10.0 ±0.8 


U.UZU_Q 






B+ 




6.9 ±2.3 


0.04 ±0.29 






B'- 


K* + TT- 


10.3 ± 1.1 


-0.25 ±0.11 






B^ 


^ K*%' 


2.4 ±0.7 


-0.15 ±0.12 




(11) 


B+ 
B+ 


^p+K'> 


q 01+0.48 

'-'•°-'--0.46 

8.0l}i 


-0.12 ±0.17 






5° 


p-K+ 


8.6+?-? 


0.15 ±0.06 








^p'K' 


rr 4+0.9 


0.01 ±0.20 


0.63l°il 



Table 1: Branching ratios, direct CP asymmetries Aqp, and mixing-induced CP 
asymmetry Scp (if applicable and measured) in S ^ K*Tr and B pK decays. 
Quoted values are taken from Ref. [8] that weight-averages individual measurements 
P[ini[IIll2l[T3l[TlEl[l6l[I3[l8l[l9l[20l[2B 

In both Fitl and Fit2, one of the theoretical parameters is the weak phase 7. We 
have two options in treating it. One possibility is to extract 7 from the B K*n 
and B —>■ pK measurements. This value can then be compared with that obtained 
from independent measurements [23] : 

7 = (66.8^^:^)° . (5) 

If the two values of 7 disagree, this will be a hint of NP. Alternatively, one can 
impose the value of 7 of Eq. ([5]) by adding a constraint to the fits. We perform 
fits using both of these options. Fitl-a and Fit2-a extract 7 from the data; Fitl-b 
and Fit2-b have the additional 7 constraint (in which the asymmetrical errors are 
averaged) . 

In all fits, there will be a hint of NP if the extracted values of the diagrams 
disagree with the SM calculations. In the QCD factorization (QCDF) approach, the 
default parameters have the following ratios of amplitudes [23] : 

\Cy/T'p\ = 0.16 ± 0.11 , and \Tp/Pi^p\ = 0.59 ± 0.12 ; 

\C'p/T{,\ = 0.20 ± 0.13 , and \T{,/Pl^y\ = 0.55 ± 0.26 . (6) 

We begin with Fitl-a. Here, the diagram P^^ is neglected. There are therefore a 
total of 11 hadronic theoretical parameters: 6 amplitude magnitudes, and 5 relative 
strong phases. In addition, we have the unknown weak phase 7. Finally, the mixing- 
induced CP asymmetry Scp{p^K^) depends also on the weak phase [3. For this, we 
impose the additional constraint (3 = (21.66lo'g7)° which is obtained mainly 
from the measurement of CP violation in B^{t) —>■ J/'^Kg and other h — > ccs decays. 
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We average the asymmetrical errors on (3. The degree of freedom {d.o.f.) is therefore 
3 in this fit. 

Before presenting the results of this fit, we note the following. There are quite 
a few solutions, many of which suggest the presence of new physics (NP). While we 
cannot be absolutely sure that one of these is not the true solution, we make the 
assumption that the NP, if present, is not enormous. If it were, it probably would 
already have been seen elsewhere. The large-NP scenarios can be eliminated by 
demanding that an acceptable solution respect certain constraints on the ratios of 
magnitudes of diagrams. Now, as we will see, the errors in the fits of the theoretical 
parameters are quite large. If we consider the full la range of these parameters, 
all potential solutions will respect the constraints. For this reason, we concentrate 
only on the central values. In particular, we require that the central values of any 
solution do not violate the ratios of Eq. ([6]) by a large amount. We impose the 
conservative constraints 

\C'ym, \C'p/T^\ < 1 , and \T'p/Pi^pl K/Kyl < 2 ■ (7) 

If any one of these is violated, the solution is excluded. There are a great many 
possible solutions. But if we restrict Xmin be reasonably small, there are only 14. 
The above constraints eliminate 12 of these. 

As we will see, one of the conclusions of our analysis is that it is important to 
reduce the experimental errors on the measurements. When this is done, the 
distribution will be deeper, and there will be far fewer minima. 

The two solutions which are retained are shown in Table [2l Both of them have 
very good fit qualities - 81% and 78%, respectively. (The quality of fit depends on 
Xmin S'lid d.o.f. individually. A value of 50% or more is a good fit; fits which are 
substantially less than 50% are poorer.) 

We now turn to Fitl-b. In this fit, the additional constraint on 7 [Eq. ([5])] is 
imposed, increasing the d.o.f. by one to 4. The assumption that any NP is not 
enormous implies again that the conditions of Eq. ([7]) are required, leading to the 
elimination of 10 of 13 possible solutions. The three that are retained are given in 
Table [31 They have good fit qualities - 64%, 58% and 56%, respectively. 

In Fit2-a, we keep the diagrams P^^ p and P^^ y. Their addition increases the 
number of hadronic theoretical parameters by 4 (2 magnitudes, 2 strong phases). 
However, the direct CP asymmetries in B ^ K*^7t~^ and B — p^K^ are no longer 
zero, leading to an increase in the number of observables by 2. The net effect is 
that the d.o.f. is reduced by 2 with respect to that of Fitl-a, i.e., it is 1. In the 
SM, it is expected that iPic pI PL p\^ Wacv I ^Lvl ~ (i'^ ia.ct, it has been argued 
in Ref. [25] that the \Pic/ Ptc\ ratio is smaller). In order to avoid too-large NP, in 
addition to the constraints of Eq. ([7]), we also require that 

\Puc,pI Ptc,p\i \Pucy I Ptcy \ < 1 • (8) 
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p' 1 

^tc,P\ 


\Pic.v\ 






0.98 


32.9ii:g 


29.0+^:^ 


13.3+^6' 


11-7+^5 


1.08 


32.4 ± 1.6 


29.3li^ 


13.61^^7^ 


12.8t«:^ 


\C'p\ 


\C'v\ 




5 pi 




1.61^:^ 


9 9+y.7 

^■^-1.6 




(128+-)° 


(-361-)° 


0.91^:^ 


2 q+a-« 

^•^-2.3 


(-161^-)° 


(-1211-)° 


(-1361-)° 










7 




(-27 ± 180)° 


(-1191-)° 


(2i.67ii:^)° 


(82.91?^-^)° 




(76 ± 180)° 




(21.65 ±0.91)° 


(91.51^^)° 



Table 2: Solutions of Fitl-a {d.o.f. = 3) to the decay observables in 5 — > i^*7r and 
B — > piiT decays. The amplitude magnitudes are quoted in units of eV. 





1 P' 1 


1 P' 1 




\T{.\ 




1.66 


32.0l{:^ 


31.711-^ 


8.61-^ 


le.ol^i 




1.99 


32.91^;^ 


27.91^:^ 


^' •'-'-16.9 


6.9lt:^ 


2.06 


33 • 9 _2 4 


97 7+2.3 
^' • '-2.6 


15.61^:^,^ 








n\ 




Sp' 


5t' 


'^•^-2.4 




(-1431-)° 


(-1381^^)° 


(-1491-)° 


2.01^:^ 




(-8l^o)° 


(1091-)° 


(-3611^)° 


4.61^:^ 


3.21^:^ 


(-15± 13)° 


(1411iD° 


(-211-)° 






Sc' 
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(791-')° 


(-1631^)° 


(21.67 ±0.91)° 


(67.2 ±4.6)° 






(-57 ± 180)° 


(-170+if)° 


(21.67 ±0.91)° 


(67.31-^)° 






(-22 ± 180)° 


(-631-)° 


(21.75 ±0.91)° 


(68.4 ±4.3)° 



Table 3: Solutions of Fitl-b (d.o.f. = 4) to the decay observables in i? — > i^*7r and 
B —>■ pK decays. The amplitude magnitudes are quoted in units of eV. 
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The imposition of all these conditions leads to the elimination of 17 of 20 possible 
solutions. The three that are kept are given in Table HI The qualities of fit are good 

- 64%, 62% and 61%, respectively. 

Finally, we have Fit2-b. This is just like Fit2-a, except that the value of 7 is 
fixed by Eq. This increases the d.o.f. by one to 2. The constraints of Eqs. 
and eliminate 8 of 10 possible solutions. The two that are retained are given in 
Table [51 The fit qualities are 88% and 67%, respectively. 

We can now examine the various fit solutions of Tables [2ll5l We separate them 
into three categories, depending on the level of disagreement with independent mea- 
surements [Eq. (^] or the predictions of the SM [Eq. (^]. There are SM fits (dis- 
crepancies at the level of < 2a), marginal fits (discrepancies at the level of ~ 3a), 
and fits which indicate NP (discrepancies at the level of ~ 5a). If we just concen- 
trate on the central values of the theoretical parameters found in the fits, we have 
all three types, shown in Table [61 Of the ten solutions in the Tables, only two are in 
agreement with the SM. In fact, six show clear signs of NP. If Table [61 represented 
the actual situation, we would conclude that B K*it and B pK decays might 
be showing signs of new physics, just like B irK. 

Unfortunately, this is not the case. In the fits, the theoretical parameters have 
errors associated with them, shown in Tables [2}[5l The errors are so large that, when 
they are taken into account, all solutions agree with the SM. This can be seen clearly 
in the following figures. 

For Fitl-a, Fig. [H shows the correlation between 7 and \C'p/Ty\ (left plot) or 
\T'p/Pp\ (right plot). In both plots, the shaded region is given by the la ranges 
of Eqs. ([5|) and IQ. For the central values of the solutions in Table [2l ((red) dots), 
we see that 7 is only marginally consistent with one of them, and inconsistent with 
the other, indicating the presence of NP. This is what Table M also states. (There 
are no inconsistencies with the ratios of diagrams.) However, these plots also show 
the allowed regions when the fit errors on the theoretical parameters are taken into 
account: the contours of Ax^ = 1 and 2 from the best solution {xLin = 0.98) are 
indicated by the solid (blue) and dashed (green) curves. We see quite clearly that 
the shaded area is mostly contained in the "solid region," and is entirely in the 
"dashed region." That is, when we include the fit errors, there is no inconsistency 

- both solutions are in agreement with the SM. 

The situation is similar in Fig. [2l Here we see the correlation between the ratios 
\Cp/Ty\ and \Cy/Tp\ for Fitl-b. As indicated in Table [6l the shaded region is 
consistent with two central values, and inconsistent with the third. However, when 
fit errors are included, the shaded area is included in the allowed regions (essentially 
completely for the "dashed region"). Once again, when the fit errors are taken into 
account, all solutions are in agreement with the SM. 

The situation is even more striking in Fig. [3l which shows the correlation between 
the weak phase 7 and the ratio \Cp/Ty \ (left plot) or the ratios \Cp/Ty \ and \Cy/Tp\ 
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P' 
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T' 
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^O.0_ig 3 
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1 ^UC^P 1 


1 P' 1 




5p' 




20.4-- 


i-^- 1 -10.7 


q q+19.1 
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15.7^^-^ 

— OU.O 


12 6+^^-^ 

lZ.O_^2.6 






6.4l^V 




10 1+^^-^ 


(-176 ±22)° 
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St' 












Suc,p\ 


(-106^-)° 


(-151+-)° 


(-33 ± 180)° 


(179 ± 180)° 


(-169 ± 180)° 




(96+i-)° 


(-69 ±180)° 


(116 ± 180)° 


(201-)° 


(-156^?^)° 


(1461-)° 


(-128 ± 180)° 


(175 ± 180)° 


(-157 ± 180)° 
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(-59 ± 180)° 


(21.66 ±0.91)° 


(73.81--)° 












(-75 ± 180)° 


(21.67 ±0.91)° 


(62.11-:^)° 












(-136 ± 180)° 


(21.67 ±0.92)° 


(116.21-:^)° 













Table 4: Solutions of Fit2-a {d.o.f. = 1) to the decay observables in 5 — > K*tt and 
S — > pK decays. The amphtude magnitudes are quoted in units of eV. 
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29.117:8 




16.7t[ti 
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P' 


5 pi 


13.1^^^;} 


q c+^i.a 
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-'-^•-'--14+ 


-'--'-•"-11.6 


(10 ±33)° 


i.2ii:- 


-L -^■--11.2 


q 1+20, 


2 






Sp' 






Sn' 






Sn' 


\Kc,p\ 


(1121-)° 


(931-^)° 


(-71 ± 180)° 


(105 ± 180)° 


(17 ± 180)° 


(1081?^)° 


(-50li-)° 


(112 ± 180)° 


(-159 ±180)° 


(-14 ± 180)° 


\Suc,v\ 


/? 
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(-77 ± 180)° 


(21.67±0.91)° 


(66.71-)° 








(115 ± 180)° 


(21.68 ±0.92)° 


(66.91^^)° 









Table 5: Solutions of Fit2-b {d.o.f. = 2) to the decay observables in S — > K*n and 
B — > pK decays. The amplitude magnitudes are quoted in units of eV. 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 o.O 0.2 0.4 0.6 0.8 1.0 

|C',/TVI |T,/PV| 

Figure 1: Fitl-a: correlation between the weak phase 7 and the ratio \C'p/Ty\ (left 
plot) or \Tp/Pp\ (right plot). The (red) dots correspond to the solutions listed in 
Table [21 The contours of Ax^ = 1 and 2 from the best solution (Xmm = 0.98) are 
indicated by the solid (blue) and dashed (green) curves. The shaded region is given 
by the Icr ranges of Eqs. ([5]) and ([6]) for comparison. 




0.0 0.2 04 0.6 0.8 1.0 

|C\./T',| 



Figure 2: Fitl-b: correlation between the ratios \C'p/Ty\ and \Cy/T'p\. The (red) 
dots correspond to the solutions listed in Table [31 The contours of Ax^ = 1 and 2 
from the best solution (Xmin — 1-66) are indicated by the solid (blue) and dashed 
(green) curves. The shaded region is given by the la ranges of Eqs. ^ and ^ for 
comparison. 
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Type of fit Table Xmm Problem 



SM 


I3l 


1.99 






m 


2.06 




marginal 


m 


0.98 






m 


0.81 


nm 


NP 


m 


1.08 






m 


1.66 


\C'v/T'p\ 




m 


0.22 


\c'yirp\ 




m 


0.24 


\C'plT'y\ 




m 


0.26 






m 


0.26 


\C'plT'y\ 



Table 6: Type of fits in Tables [2ll5l using only the central values of the theoretical 
parameters. 

(right plot) for Fit2-a. Taking the two plots together, one sees that all three central 
values ((red) dots) are inconsistent with the shaded region, that is, they all indicate 
NP. However, the shaded region is completely inside the "solid region." Thus, the 
inclusion of the fit errors removes all inconsistencies, so that all solutions are in 
agreement with the SM. 

Finally, Fig. H] shows the correlation between the ratios |Cp/Ty| and \CylT'p\ for 
Fit2-b. Despite the fact that the central values show some discrepancies with the 
shaded region, the inclusion of the fit errors removes them, so that all solutions are 
in agreement with the SM. 

From the figures, one sees that, when one includes the errors from the fits on the 
theoretical parameters, all solutions are consistent with the SM, even if the central 
values are not. The reason for this is that the fit errors are quite large. And this 
is related to the fact that the experimental errors on the B — >■ K*!^ and B —>■ pK 
observables are also large, especially those related to the CP asymmetries. These 
errors are considerably larger than those in B ^ Ktt. At present, it is not possible to 
clearly see if NP is present in these 6 — > s VP decays. Hopefully, in the future, it will 
be possible to increase the precision of the B K*tt and B pK measurements, 
so that we can further test the irK puzzle. 

In summary, we have performed fits to -B K*7c and B pK decays within 
the standard model (SM). In particular, we have employed the relations between 
electroweak-penguin and tree amplitudes to couple together the two types of decay 
modes, which are otherwise separate. Our fits are of two types: those in which the 
weak phase 7 is extracted from the i?-decay data, and those in which an independent 
determination of 7 [Eq. ([5])] is imposed. We begin our analysis by neglecting P^^. We 
find that the central value of 7 tends to be larger when no independent determination 
is imposed (Fitl-a). When Eq. ^ is taken into account, two of the solutions are 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.2 0.4 0.6 0.8 1.0 

|C',/TVI ICV/TVI 



Figure 3: Fit2-a: correlation between the weak phase 7 and the ratio \C'p/Ty\ (left 
plot) or the ratios \C'p/Ty\ and \C'y/T'p\ (right plot). The (red) dots correspond 
to the solutions listed in Table HI The contour of Ax^ = 1 from the best solution 
ixlnin — 0.22) is indicated by the solid (blue) curve. The shaded region is given by 
the la ranges of Eqs. ([5]) and ([6]) for comparison. 

consistent with the SM expectations, while the other has a somewhat large \Cy/T'p\ 
ratio (Fitl-b). Next we include P^^ in the fits. Without using the external constraint 
on 7, all three solutions have a problem with too-large \C' /T'\ ratios. In addition, 
one of these gives a too-large value for 7 (Fit2-a). They could indicate the presence 
of new physics (NP). After imposing Eq. ([5]), there are two solutions. One of these 
still has the problem of a large |C"/T'|, while the other is improved (Fit2-b). 

These problems with the theory parameters potentially show some hints of NP. 
However, the fits also yield large errors on these parameters. There is no serious 
discrepancy with the predictions of the SM when these errors are taken into account. 
In order to make more decisive conclusions about NP m. B ^ K*7i and B pK, the 
uncertainties on the measurements of observables have to be reduced, particularly 
for the CP asymmetries. Such efforts will be paid off by shedding additional light 
on the B nK puzzle. 
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Taiwan (CC) and by NSERC of Canada (DL). 



10 




0.8 1.0 



Figure 4: Fit2-b: correlation between the ratios \C'plTy\ and \C'y/T'p\. The (red) 
dots correspond to the solutions listed in Table [51 The contours of Ax^ = 1 and 2 
from the best solution (Xmm — 0.26) are indicated by the solid (blue) and dashed 
(green) curves. The shaded region is given by the la ranges of Eqs. ^ and ^ for 
comparison. 
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